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Abstract 

In this article we will present pure two-dimensional analytic solutions for the coupled non- 
compressible Newtoniain Navier-Stokes — with Boussinesq approximation — and the heat con¬ 
duction equation. The system was investigated from E.N. Lorenz half a century ago with Fourier 
series and pioneered the way to the paradigm of chaos. We present a novel analysis of the same 
system where the key idea is the two-dimensional generalization of the well-known self-similar 
Ansatz of Barenblatt which will be interpreted in a geometrical way. The results, the pressure, 
temperature and velocity fields are all analytic and can be expressed with the help of the error 
functions. The temperature field has a strongly damped oscillating behavior which is an interesting 
feature. 
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I. INTRODUCTION 


The investigation of the dynamics of viscous fluids has a long past. Enormous scientific 
literature is available from the last two centuries for fluid motion even without any kind 
of heat exchange. Thanks to new exotic materials like nanotubes, heat conduction in solid 
bulk phase (without any kind of material transport) is an other quickly growing independent 
research area as well. The combination of both processes are even more complex which 
lacks general existence theorems for unique solutions. The most simple way to couple these 
two phenomena together is the Boussinesq [l| approximation which is used in the held of 
buoyancy-driven how (also known as natural convection). It states that density diherences 
are sufficiently small to be neglected, except where they appear in terms multiplied by g, 
the acceleration due to gravity. The main idea of the Boussinesq approximation is that the 
diherence in inertia is negligible but gravity is sufficiently strong to make the specihc weight 
appreciably diherent between the two huids. When the Boussinesq approximation is used 
than no sound wave can be described in the huid, because sound waves move via density 
variation. 

Boussinesq hows are quite common in nature (such as oceanic circulations, atmospheric 
fronts or katabatic winds), industry (fume cupboard ventilation or dense gas dispersion), 
and the built environment (like central heating, natural ventilation). The approximation is 
extremely accurate for such hows, and makes the mathematics and physics much simpler 
and transparent. 

The advantage of the approximation arises because when investigation a how of, say, 
warm and cold waters with densities pi and p 2 are considered, the diherence Ap = pi — P 2 is 
negligible and one needs only a single density p. It can be shown with the help of dimensional 
analysis, under these circumstances, the only sensible way that acceleration due to gravity 
g should enter into the equations of motion is in the reduced gravity g' = g{pi — p 2 )- 
The corresponding dimensionless numbers for such hows are the Richardson and Rayleigh 
numbers. The used mathematics is therefore much simpler because the density ratio (pi/p 2 
a dimensionless number) is exactly one and does not ahect the features of the investigated 
how system. 

In the following we analyze the dynamics of a two-dimensional viscous huid with addi¬ 
tional heat conduction mechanism. Such systems were hrst investigated by Boussinesq [l| 
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and Oberbeck [2j in the nineteenth century. Oberbeck used a hnite series expansion. He 
developed a model to study the heat convection in fluids taking into account the flow of the 
fluid as a result of temperature difference. He applied the model to the normal atmosphere. 
More than half a century later Saltzman 3| tried to solve the same model with the help 


of Fourier series. At the same time Lorenz 


^ analyzed the solutions with computers and 


published the plot of a strange attractor which was a pioneering results and the advent of 
the studies of chaotic dynamical systems. The literature of chaotic dynamics is enormous 


fl. 


but a modern basic introduction can be found in _ 

Later till to the first beginning years of the millennium jd] Lorenz analyzed the final first 
order chaotic ordinary differential equation(ODE) system with different numerical methods. 
This ODE system becomes an emblematic object of chaotic systems and attracts much 
interest till today jb]. 

On the other side critical studies came to light which go beyond the simplest truncated 
Fourier series. Curry for example gives a transparent proof that the finite dimensional 
approximations have bounded solutions [^. Musielak ef a/ js] in three papers analyzed large 
number of truncated systems with different kinds and found chaotic and periodic solutions 
as well. The messages of these studies will be shortly mentioned later. 

In our study we apply a completely different investigation approach, namely the two- 
dimensional generalization of the self-similar Ansatz which is well-known for one dimension 


from more than half a century jb -If l| . This generalized Ansatz was successfully app 


three dimensional compressible and non-compressible Navier-Stokes equations 12 


ied to the 
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us in the last years. We investigated one dimensional Euler equations with heat conduction 
as well which can be understood as the precursor of the recent study. 

To our knowledge this kind of investigation method was not yet applied to the Oberbeck- 
Boussinesq (OB) system. In the next section we outline our theoretical investigation together 
with the results. The paper ends with a short summary. 


II. THEORY AND RESULTS 

We consider the original partial differential equation(PDE) system of Saltzman to 
describe heat conduction in a two dimensional viscous incompressible fluid. In Cartesian 
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coordinates and Eulerian description these equations have the following form: 



0 , 


0 , 


0 , 


du dw 
dx dz 


0 , 


(1) 


where u,w, denote respectively the x and z velocity coordinates, Ti is the temperature 
difference relative to the average (Ti = T — Tav) and P is the scaled pressure over the 
density . The free physical parameters are z/, e, G, k kinematic viscosity, coefficient of volume 
expansion, acceleration of gravitation and coefficient of thermal diffusivity. (To avoid further 
misunderstanding we use G for gravitation acceleration and g which is reserved for a self¬ 
similar solution.) The hrst two equations are the Navier-Stokes equations, the third one is the 
heat conduction equation and the last one is the continuity equation all are for two spatial 
dimensions. The Boussinesq approximation means the way how the heat conduction is 
coupled to the second NS equation. Chandrasekhar jl^ presented a wide-ranging discussion 
of the physics and mathematics of Rayleigh-Benard convection along with many historical 
references. 

Every two dimensional flow problem can be reformulated with the help of the stream 
function T via m = and v = —which automatically fulhlls the continuity equation. 
The subscripts mean partial derivations. After introducing dimensionless quantities the 
system of ([T]) is reduced to the next two PDEs 


+ '^yy)t + ^x{:^xxz T yyz) - + ^zzx) - 

(y{Qx - ^xxxx - ^zzzz - "i^xxzz) = 0, 

Qt + "^x^z — — {0XX + (^zz) = 0 , 


xxxx 


( 2 ) 


where 0 is the scaled temperature, a = v /kis the Prandtl Number and R = is the 

Rajdeigh number and H is the height of the fluid. A detailed derivation of ([2]) can be found 

in yl. 


All the mentioned studies in the introduction, investigated these two PDEs with the 
help of some truncated Fourier series, different kind of truncations are available which result 
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different ordinary differential equation(ODE) systems. The derivation of the final non-linear 
ODE system from the PDE system eatt be fouad m the otigiaa. papers BQ. Betge et a/, 
[Iq contains a slightly different development of the Lorenz model equatios, and in addition, 
provides more details on how the dynamics evolve as the reduced Rayleigh number changes. 
The book of Sparrow gives a detailed treatment of the Lorenz model and its behavior 
as well. Hilborn presents the idea of the derivation in a transparent and easy way. 
Therefore, we do not mention this derivation in our manuscript. 

Some truncations violates energy conservation Q| and some not. Roy and Musiliak js] in 
his exhausting three papers present various energy-conserving truncations. Some of them 
contain horizontal modes, some of them contain vertical modes and some of them both kind 
of modes in the truncations. All these models show different features some of them are 
chaotic and some of them - in well-defined parameter regimes - show periodic orbits in the 
projections of the phase space. This is a true indication of the complex nature of the original 
flow problem. It is also clear that the Fourier expansion method which is a two hundred 
year old routine tool for linear PDFs fails for a relevant non-linear PDF system. 

Therefore, we apply another investigation method which is common for non-linear PDFs. 
At hrst we introduce the two dimensional generalization of the self-similar Ansatz 


v{x,t)=t-fl^^y.= t-f{v) (3) 

where v{x,t) can be an arbitrary variable of a PDF and t means time and x means spatial 
dependence. The similarity exponents a and (3 are of primary physical importance since 
a represents the rate of decay of the magnitude v{x,t), while (3 is the rate of spread (or 
contraction if < 0 ) of the space distribution for t > 0. The most powerful result of this 
Ansatz is the fundamental or Gaussian solution of the Fourier heat conduction equation 
(or for Pick’s diffusion equation) with a = [3 = 1/2. These solutions are exhibited on 
Figure 1. for time-points ti < ^ 2 - This transformation is based on the assumption that 
a self-similar solution exists, i.e., every physical parameter preserves its shape during the 
expansion. Self-similar solutions usually describe the asymptotic behavior of an unbounded 
or a far-held problem; the time t and the space coordinate x appear only in the combination of 
f{x/t^). It means that the existence of self-similar variables implies the lack of characteristic 
lengths and times. These solutions are usually not unique and do not take into account the 
initial stage of the physical expansion process. It is also transparent from (|3]) that to avoid 
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singularity at t = 0 the following transformation t = t + to is valid. 

There is a reasonable generalization of ([3]) in the form of v{x,t) = h(t) ■ f[x/g(t)], where 
h{t),g{t) are continuous functions. The choice of h{t) = g{t) = y'fo ~ t is called the blow¬ 
up solution, which means that the solution becomes inhnity after a well-dehned hnite time 
duration. 

These kind of solutions describe the intermediate asymptotic of a problem: they hold 
when the precise initial conditions are no longer important, but before the system has 
reached its hnal steady state. For some systems it can be shown that the self-similar solution 
fulhlls the source type (Dirac delta) initial condition. They are much simpler than the full 
solutions and so easier to understand and study in different regions of parameter space. A 
hnal reason for studying them is that they are solutions of a system of ODEs and hence do 
not suffer the extra inherent numerical problems of the full PDEs. In some cases self-similar 
solutions helps to understand diffusion-like properties or the existence of compact supports 
of the solution. 

Let’s introduce the two dimensional generalization of the self-similar Ansatz ([3]) which 
might have the general form of 


v{x,z,t) = t "/ 


F{x, z] 


( 4 ) 


where F{x, z) could be understood as an implicit parametrization of a one-dimensional space 
curve with continuous hrst and second derivatives. In our former studies [l2, 13] we explain 
in heavy details that for the Navier-Stokes type of non-linearity unfortunately only the 
F{x, z) = X + z + c function is valid in Cartesian coordinates which is a straight line. It 
basically comes from the symmetry properties of the left-hand of the NS equation. Only 
this function fulhlls the following relation = u^- (Other locally orthogonal coordinate 
systems e.q. spherical are not investigated yet.) 

We may investigate both dynamical systems, the original hydrodynamical ([I]) or the other 


one (ED which is valid for the stream 


’unctions. 


Similar to the former studies 


, try to solve the PDEs for the dimensionless stream 


and temperature functions in the form of 


4^ = t "/(?/), 0 t 'h(r/), 7] = (5) 

Unfortunately, after some algebra it becomes clear that the constraints which should hx 
the values of the exponents become contradictory, therefore no unambiguous ODE can be 
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FIG. 1: A self-similar solution of Eq. (I3|) for ti < t 2 - The presented curves are Gaussians for 
regular heat conduction. 

derived. This means that the PDE of the stream function and the dimensionless temperature 
do not have self-similar solutions. In other words these functions have no such a diffusive 
property which could be investigated with the self-similar Ansats, which is a very instructive 
example of the applicability of the trial function of (|5]). Our experience shows that, most of 
the investigated PDEs have a self-similar ODE system and this is a remarkable exception. 

Now investigate the original hydrodynamical system with the next Ansatz 

u{ri) = f“/(r 7 ), w{ri) = t-^g{ri), P{v) = t~"h{ri), Ti{ri) = (6) 

where the new variable is rj = {x + z)/t^. All the hve exponents a,P,6,e,u are real 
numbers. (Solutions with integer exponents are the self-similar solutions of the hrst kind 
and sometimes can be obtained from dimensional considerations.) The f,g,h,l objects are 
called the shape functions of the corresponding dynamical variables. 

After some algebraic manipulations the following constrains are hxed among the self¬ 
similarity exponents : a = 6 = f3 = 1/2, e = 1 and u = 3/2 which are called the universality 
relations. At this point it is worth to mention that now all the exponents have a £x numerical 
value which simplihes the structure of the solutions. There is no free exponential parameter 
in the original dynamical system, like an exponent in the equation of state. As an example 
we mention one of our former study where the compressible NS equation was investigated 
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FIG. 2: Different shape functions of the temperature Eq. ([9]) as a function of r] for different thermal 
diffusivity. The integration constants are ci = C 2 = 1 the same for all the three curves. The solid 
the dashed and the dotted lines are for k = 1,2,5, respectively. 

with a free parameter which described different materials. 

These universality relations dictate the corresponding coupled ODE system which has 
the following form of 


- f - ^ + //' + 9/' + fc' - ii'f" = 0, 

~2 ~ ~ ~ 2z/5f" = 0, 

-f - Y + /'' + 9i' - 2ri" = 0, 

/' + 9' = 0. (7) 


From the last (continuity) equation we automatically get the f + g = c and /" + g" = 0 
conditions which are necessary in the following. 

Going through a straightforward derivation the next single ODE for the shape function 
of the temperature distribution can be separated 

I'n V 
~2 ~2 ^ 


8 







250 - 


200 - 


150 - 

h(ri) 


100-1 



0 10 20 30 




FIG. 3: Different shape functions of the pressure Eq. (12) as a function of r] for different thermal 
diffusivity. The integration constants are taken ci = C 2 = 1 for all the three curves. We fixed the 
value of eG = 1 as well. The solid, the dashed and the dotted lines are for k = 1,2,5 numerical 
values, respectively. 


The solution is 


I = Cl 


4erfi 





^ + 4\/k?7 


-k C2e 8k (4/t — 


(9) 


where ci,C 2 are free integration constants. The erh means the imaginary error function 
dehned via the integral 2/y^ exp(x^)dx for more details see 19|. It is interesting, that 
the temperature distribution is separated from the other three dynamical variables an does 
not depend on the viscosity coefficients as well. We may say, that among the solution 
obtained from the self-similar Ansazt the temperature has the highest priority and this 
quantity dehnes the pressure and the velocity held. That is a remarkable feature. In a 
ormer study, where the one-dimensional Euler system was investigated with heat conduction 
3 we found the opposite property, the density and the velocity held were much simpler 
than the temperature held. Figure 2 presents different shape functions of the temperature 
for different thermal diffusivity values. The hrst message is clear, the larger the thermal 
diffusivity the larger the shape function of the temperature distribution. A detailed analysis 
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of Eq. (9) shows that for any reasonable k and c valnes the main property of the fnnction 
is not changing - has one global maximnm and minimnm with a strong decay for large rjs. 
A second remarkable feature is the single oscillation which is not a typical behavior for self¬ 


similar solutions. We investigated numerous non-linear PDE systems till today 
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14| some 


of them are even not hydrodynamical 2^ and never found such a property. This analysis 
clearly shows that at least the temperature distribution in this physical system has a single¬ 
period anharmonic oscillation. For a hxed time value and a well-chosen z the difference 
of values of r] where 1{t]) yields a minima and a maxima correspondsto that Ax at which 
the temperature (and density) fluctuation may start the Behnard convection. To go a step 
further we may calculate the Fourier transform of the shape function, /(p) Eq. dH]) to study 
the spectral distribution. (An analytic expression for the Fourier transform is available, 
which we skip now.) The hrst term (which is proportion to Ci) becomes a complex function, 
however the general overall shape remains the same, a single-period anharmonic oscillation 
with a global minimum and maximum like on Figure 2. Of course, the zero transition of the 
function depends on the value of k. The second term of the Fourier transformed function 
which is proportional to C 2 remains a Gaussian which is not interesting. 

For completeness we give the full two dimensional temperature held as follows 


Ti{x, z, t) = Cl t 

C2 t 


Aerfi 


_n/n _Oi+£r 

Snt 


X + z 

(x + zy 

Ak - 


( K 


(a;-|-z)^\ _(£±£d A^/K{x + z) 




At 


Shtt 


tV2 


t 


+ 

( 10 ) 


The shape function of the pressure held can be obtained from the temperature shape function 
via the following equation 


h' = 


eGl 


( 11 ) 


with a similar solution to 


h = Cl 


2K\/^eG ■ erfi ( I 


+ C22eGKT]e s-. + cs. 


( 12 ) 


this can be understood that the derivative of the pressure is proportional to the temperature. 
With the known numerical value of the exponent e = 1 the scaled pressure held can be 
expressed as well P{x,z,t) = t~^h{[x + z]/t~^^y. Note, the diherence between the u and 
the e exponents, which are responsible for the diherent asymptotic decays. The temperature 
held has a stronger damping for large rj than the pressure held. (It is worth to mention that 
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for the three dimensional NS equation, without any hea 


the pressure term is also different to the velocity held 


exchange the decay exponent of 


la- 


At last the ODE for the shape function of the velocity component z reads 


Aug" + g'rj + g + eGl = 0 


(13) 


which directly depends on the temperature on l{ri) and all the physical parameters u, e, G, k, 
of course. In contrast to the pressure and temperature held there is no closed solutions 
available for a general parameter set. The formal, most general solution is 


v~ _nz. I f 1 

g = C 2 e f — 


Cl 


AeGKC 2 'r)e — AeGny/^cier fi {rje ] e^’' 



where ci and C 2 are the recent integration constants. Note, that the integral can be analyt¬ 
ically evaluated if and only if u = k which is a great restriction to the physical system. We 
skip this solution now. The other way is to hx ci = 0 and let k and g, free. The solution 
has the next form of 


r I V 

g = Cie »‘'erj 


I 2 \ _ AeGc2K^e «« 

— ] 1 + C2e 8 ^- 

U K — U 


(15) 


Note, that now the u ^ k condition is obtained. The Ci and C 2 are the recent integration 
constants as above, it is interesting that if both of them are set to zero, the solution is still 
not trivial. For a physical system the kinematic viscosity i/ > 0 is always positive, therefore 
in the case of ci 7 ^ 0 the solution becomes complex. Figure 4 shows the shape function of the 
z velocity component. It is clear that the real part is a Gaussian function and the complex 
part is a Gaussian distorted with an error function, which is an interesting hnal result. 


In the literature we can End system which shows similarities like the work of Ernst 


2l| 


who presented a study where a the asymptotic normalized velocity autocorrelation function 
calculated from the linearized Navier-Stokes equation has an error function shape. 


III. SUMMARY 

We investigated the classical OB equation which is the starting point of countless dy¬ 
namical and chaotic systems. Instead of the usual Fourier truncation method we applied 
the two-dimensional generalization of the self-similar Ansatz and found a coupled non-linear 
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FIG. 4: The shape functions of the z velocity field g{ri) Eq. ()15p as a function of rj. The solid line 
is the real and the dotted is the complex part. All the integration constants are taken ci = 62 = 
C 2 = 1. The physical constant eG = 1 as well. The n = 0.04 and v = 0.8. 

ODE system which can be solved with quadrature. The main result is that even these kind 
of solutions - build up from error functions - show some oscillating behavior. The result¬ 
ing expressions show the decay of temperature and pressure fluctuations and velocity held 
in time. Due to our knowledge certain parts of the climate models are based on the OB 
equations therefore our results might be an interesting sign to climate experts. 
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